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Abstract 

In this paper, we develop an 0 ((to log fc)MSF(n, m, l))-time algorithm to find a half¬ 
integral node-capacitated multiflow of the maximum total flow-value in a network with n 
nodes, m edges, and k terminals, where MSF(n', m', 7 ) denotes the time complexity of solv¬ 
ing the maximum submodular flow problem in a network with n' edges, m! edges, and the 
complexity 7 of computing the exchange capacity of the submodular function describing 
the problem. By using Fujishige-Zhang algorithm for submodular flow, we can find a max¬ 
imum half-integral multiflow in 0{mn^ log k) time. This is the first combinatorial strongly 
polynomial time algorithm for this problem. Our algorithm is designed on the basis of a 
developing theory of discrete convex functions on certain graph structures. Applications 
include “ellipsoid-free” combinatorial implementations of a 2 -approximation algorithm for 
the minimum node-multiway cut problem by Garg, Vazirani, and Yannakakis. 

Keywords: Node-capacitated multiflow, discrete convex analysis, submodular flow, node¬ 
multiway cut 


1 Introduction 

A node-capacitated undirected network is a quadruple N = {V,E,S,c) of node set V, (undi¬ 
rected) edge set E, a specified subset S of nodes, called terminals, and a nonnegative integer¬ 
valued node capacity c : V \ S ^ Z+ on nonterminal nodes. An S-path is a path connecting 
distinct terminals. A (node-capacitated) multiflow is a pair {V, A) of a set V of S'-paths and 
a flow-value function A : "P —)• R+ satisfying the node-capacity constraint: 

X KP)<c{i) {iev\s). ( 1 . 1 ) 

Per: ieV(P) 

The total flow-value of a multiflow / = (P, A) is defined as YlpeV ^ multiflow is called 

maximum if it has the maximum total flow-value among all possible multiflows. A multiflow 
/ = (P, A) is said to be integral if A is integer-valued, and half-integral if 2A is integer-valued. 

In this paper, we address the problem of finding a maximum multiflow in a node-capacitated 
network. This multiflow problem appeared in the work by Garg, Vazirani, and Yannakakis |T2] 
on an approximation algorithm for node-multiway cut. In fact, the LP-dual of our multiflow 
problem is a natural LP-relaxation of the minimum node-multiway cut problem; see also m 
Section 19.3]. They showed that this LP-dual always has a half-integral optimal solution. The 
half-integrality of the primal problem, i.e., the existence of a half-integral maximum multiflow. 
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was later shown by Pap |24l [25] . He also showed that a half-integral maximum multiflow can 
be found in strongly polynomial time. 

In these works, the polynomial time solvability depends on the use of the ellipsoid method. 
Thus it is natural to seek a combinatorial polynomial time algorithm. For the case of unit 
node-capacity (c(i) = 1 for all i G V\S), Babenko [2] developed a combinatorial 0{mv?) time 
algorithm to find a half-integral maximum multiflow, where n is the number of nodes and m 
is the number of edges. For general node-capacity, Babenko and Karzanov [3] developed a 
combinatorial weakly polynomial time algorithm to find a half-integral maximum multiflow. 
Their algorithm runs in 0{MF{n, m, C)n^ \og^ nlogC) time, where MF(n, 771,(7) is the time 
complexity of solving the max-flow problem in a network with n nodes, m edges, and the 
maximum edge-capacity C. 

The main result of this paper is the first combinatorial strongly polynomial time algorithm 
to solve the maximum node-capacitated multiflow problem. Our algorithm uses, as a sub¬ 
routine, an algorithm of solving the maximum submodular flow problem; see [9l Section 5.5 
(c)]. Let MSF( 77 , 777 , 7 ) denote the time complexity of solving the maximum submodular flow 
problem on a network with n nodes, m edges, and the time complexity 7 of computing the 
exchange capacity of the submodular function describing the problem. 

Theorem 1.1. There exists an 0{{mlogk)M.SF(n,m, l))-time algorithm to find a half-integral 
maximum multiflow and a half-integral optimal dual solution in a network ofn nodes, m edges, 
and k terminals. 

The current fastest maximum submodular flow algorithm is the push-relabel algorithm due 
to Fujishige and Zhang |llj of the time complexity 0 ( 77 ^ 7 ); see the survey [10] on submodular 
flow algorithms. Thus we can solve the problem in 0{mn^ log A;) time. 

Application 1: Node-multiway cut. A node-multiway cut is a subset X C H \ S' of 
nonterminal nodes such that the deletion of X makes every pair of distinct terminals un¬ 
reachable, or equivalently, X meets every S-path. The capacity of a node-multiway cut X is 
defined as Yliex '^(^)- minimum node-multiway cut problem asks to find a node-multiway 
cut with the minimum capacity. This well-known NP-hard problem is naturally formulated 
as the following { 0 , l}-integer program: 

Minimize c{i)w{i) 

iev\s 

subject to w : V \ S ^ {0,1}, 

w{i) > 1 (every S-path P). (1.2) 

ieViP)\s 

The natural LP-relaxation obtained by relaxing tc : P \ S —>■ {0,1} into w : V \ S ^ R+ 
is nothing but the LP-dual of our multiflow problem. As mentioned above, Garg, Vazirani, 
and Yannakakis [T2| proved that a half-integral optimal LP solution w* : P \ S —)• {0,1/2,1} 
always exists, and is obtained from any optimal LP solution by a simple rounding procedure; 
see m Section 19.3]. Then the set of nodes i with w*{i) > 1/2 is a 2-approximation solution 
of the minimum node-multiway cut problem. This rounding algorithm needs an optimal LP 
solution, which is now obtained by our algorithm. To the best of our knowledge, this is the first 
combinatorial strongly polynomial time implementation of the 2 -approximation algorithm. 

Just recently, Chakuri and Madan |5] devised a simple method to round any feasible LP 
solution into a multiway cut of capacity within factor 2. Combining this rounding method with 
a fast FPTAS for multiflow (e.g., [I3]), they obtain a considerably fast (2 -|- e)-approximation 
algorithm (with running time dependent of 1 /e). 
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Application 2: Integral multiflow. Our algorithm is also useful in the problem of finding 
a maximum integral multiflow. This problem is a capacitated version of openly-disjoint S- 
paths packing problem considered by Mader [2T] (that corresponds to the case of c{i) = 1 
{i G V\S)). Pap [241125] established the strongly polynomial time solvability of the maximum 
integral multiflow problem. The first step of his algorithm is to find a maximum half-integral 
multiflow. The second step is to construct and solve an instance of the openly-disjoint S-path 
packing problem (on the graph with size polynomial in the numbers of edges and nodes in 
the original network). Finally, combining the integer part of the half-integral multiflow with 
a solution of the packing problem, one obtains a maximum integer multiflow. The second 
step can be done by several combinatorial polynomial time algorithms, including [6] and [261 
Section 73.1a]. Our algorithm can be used in the first step, and makes the whole algorithm 
fully combinatorial. 


Outline. Let us outline our algorithm and the ideas behind it, as well as the structure of 
the paper. Our algorithm is designed on the basis of the following two ingredients. One is a 
combinatorial duality theory for a class of node-capacitated multiflow problems [l5] . The other 
is a developing theory of discrete convex functions on certain graph structures [m [izi HU, 
which aims to extend concepts in Discrete Convex Analysis (DCA) (Murota [22]) to tackle 
further various combinatorial optimization problems beyond network flows, matroids, and 
submodular functions. We will utilize these theories in a self-contained way. 

In Section following [15] we formulate the dual of our multiflow problem as a facility 
location problem on a tree. This formulation gives a fruitful combinatorial interpretation of 
the LP-dual problem (1.2), and brings a simple combinatorial algorithm to find a half-integral 
optimal multiflow from a given dual optimum, under a certain nondegeneracy assumption. 
We will deal with a perturbed problem satisfying this nondegeneracy assumption. Our goal 
is to solve this perturbed problem efficiently. We will see that the location problem is further 
formulated as an optimization over a certain discrete structure, and the objective is an L- 
convex function on a Euclidean building in the sense of |18j . This class of discrete convex 
functions shares many analogous properties with L-convex functions in DCA. In particular, as 
in the case of DCA, there is a natural descent algorithm, called the steepest descent algorithm, 
to minimize our L-convex function g. For each point x, the steepest descent algorithm chooses 
a point y {steepest direction) from a discrete neighborhood of x with smallest g{y). If g{y) = 
g{x), then x is guaranteed to be optimal. Otherwise, i.e., g{y) < g{x), replace x by y, and 
repeat. 

In Section]^ we will implement this conceptually simple algorithm. We will prove that in 
our case a steepest direction at each point can be found by solving one maximum submodular 
flow problem. This part is the heart of our analysis. As a consequence, we obtain an algorithm 
in a simple form as follows: 


1. From a dual solution x, construct and solve an instance of the maximum submodular 
flow problem. 

2. If the minimal minimum cut consists only of the source, then x is optimal, and an 
optimal multiflow is constructed from any maximum submodular flow. Otherwise the 
minimal minimum cut gives a steepest direction y of the neighborhood at x. Replace x 
by y, and go to 1. 


Our maximum submodular flow problem is defined by a disjoint sum of submodular functions 
on 6-element sets. This enables us to compute the exchange capacity in constant time. More¬ 
over, the number of iterations is estimated by the geodesic descent property (Theorem |3.5[ ) of 
the steepest descent algorithm. This intriguing property says that a trajectory of the algo¬ 
rithm forms a geodesic to optimal solutions with respect to a certain /oo-metric on the domain. 
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We know in advance the range where an optimum exists, and the diameter of the range is 
bounded by 0(m log A:) relative to the above metric. Consequently the number of iteration is 
bounded by 0{mlogk). 

It should be noted that our algorithm design includes an interesting new technique of 
reducing bisubmodularity to submodularity. This technique and related arguments, including 
basics on submodularity, are summarized in Section Actually step 1 of the above algorithm 
is essentially the feasibility check of a bisubmodular flow problem, that is, hnding a (fractional) 
bidirected flow with the flow-boundary constrained to a bisubmodular polyhedron. This 
seemingly natural class of problems has not been well-studied so far. It is well-known that 
(fractional) bidirected flows are easily manipulated by ordinary flows in a skew-symmetric 
network obtained by doubling nodes and edges. We generalize this doubling construction to 
bisubmodular functions. We give a condition for a bisubmodular function / to be extended 
to a submodular function /' on a larger set, so that the bisubmodular polyhedron of / is a 
projection of the base polyhedron of f. We show that a certain bisubmodular function on 
a 3-element set, which represents the flow-conservation and the node-capacity constraints on 
a node of degree 3, has such a submodular extension on a 6-element set. Our bisubmodular 
flow problem is described by the disjoint sum of these bisubmodular functions, and can be 
reduced to the submodular flow problem as mentioned above. 

2 Preliminaries 

Notation. Let R, R+, Z, and Z+ denote the sets of reals, nonnegative reals, integers, and 
nonnegative integers, respectively. The infinity element oo is treated as x < oo and x-l-oo = oo 
for X G R. The set of all functions from a set V to a. set R is denoted by RY . For a function 
: F —)■ R U {oo}, let dom^r := {x G F | g{x) < oo}. For a function v G R'^ and a subset 
X C V, let v{X) denote The function value v{i) will also be denoted by vi if 

no confusion occurs. For a (directed or undirected) graph G = (F, E), an edge from i to j is 
denoted by ij. For a subset X of nodes, let 5X denote the set of all edges leaving X. For an 
undirected graph T with a specified edge-length, let d = dp denote the shortest path metric 
on the vertex set with respect to the edge-length. In the following, graphs or networks are 
supposed to have no multiple edges and loops. 

Signed set and transversal. A signed set U is the product F x {-|-, —} of a set F and the 
sign {+,—}. Elements (i,-|-) and (i, —) of U are simply denoted by Y and i~, respectively. 
The signed extension of a set F is defined as the signed set F x {+, —} and is denoted by 
F=t. For F C F, let y+ := {i+ | i G Y} and F” := ji" | i G F}. Also for U C F=t, let 
[/+ := {i+ I G 17} and U~ := {i~ \ i~ G 17}. A subset X of F^ is called transversal if 
\X n {i+, i“}| < 1 for all i G F, and is called eo-transversal if \X n {i+, i“}| > 1 for all i G F. 
For X C F^ let A denote the transversal obtained from X by deleting all {Y with 
{Y C A, and let A denote the co-transversal obtained from A by adding all {Y,i~} 
with {Y,i~} n A = 0. For u G F^, define u hy u := Y if u = i~ and u := i~ if u = Y. 

Skew-symmetric network. A skew-symmetric network (see e.g., El) is a directed network 
on a signed set such that an edge uv exists if and only if an edge vu exists, and the (lower 
and upper) capacities of edges uv and vu are the same. A skew-symmetric network is often 
useful for dealing with undirected objects. 

2.1 Submodular flow 

Here we summarize basics on submodular functions and submodular flows; see HI [g do] for 
further details. A submodular funetion on a set F is a function p defined on 2^ satisfying 
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p{X) + p{Y) > p{X n y) + p{X U y) for X, y c y. Let p be a submodular function on V with 
/?(0) = 0. The base polyhedron B{p) is defined to be the set of all vectors x G R'^ satisfying 
x{X) < p{X) for X C y and xiy) = p{V). For x G B{p) and a pair (i, j) of distinct elements 
of y, the exchange capacity k(x; i,j) at x is defined by 

Kix;i,j) := max{a G R+ | x + a{xi - Xj) £ B{p)}, 

where Xi is the i-th unit vector defined by Xi(j) := 1 if z = j and Xi(j) ■= 0 otherwise. 

We next introduce submodular flows. Let X be a directed network with vertex set V, 
edge set A, edge-capacity c : A ^ R+, and terminals s,t G V. The set of nonterminal 
nodes is denoted by U{:= V \ {s,f}). We are given a submodular function /o : 2^ —)> R with 
/?(0) = p{U) = 0. For a function 9 ? : A —)• R, let Xp denote the function on V defined by 

Xp{i) := ^^{<y 9 (e) I e G e enters 1 } — \ e G A\ e leaves i}. 

Namely Xp{i) represents the excess of p at node i. An {s,t)-flow, or simply, a flow is a 
function p : A ^ R 4 . satisfying 

0 < p{e) < c(e) (e G A), Xp\u G B{p), 

where (•)|[/ means the restriction to U. The flow-value of a flow p is defined as Xp{t){= 
—Xp{s)). An (s,t)-cut X is a subset of nodes containing s and not containing t. The cut 
capacity of X is defined as 

c(5X) + p(X\{s}). (2.1) 

An (s, f)-flow is called maximum if it has the maximum flow-value, and an (s, f)-cut is called 
minimum if it has the minimum capacity. 

Theorem 2.1 (see O Theorem 5.11]). The maximum flow-value of an {s,t)-flow is equal to 
the minimum cut-capacity of an {s,t)-cut X. If c and p are both integer-valued, then there 
exists an integer-valued maximum flow. For any maximum flow p, the set of nodes reachable 
from s in the residual network of p is the unique minimal minimum {s,t)-cut. 

Here the residual network of is a directed network on V constructed as follows: 

For each edge ij in N with p{ij) < c{ij), add an edge ij to Np. For each edge ij in N with 
0 < p{ij), add an edge ji to Np. For each pair of distinct nodes i, j in U with K{Xp\if', i, j) > 0, 
add an edge ij to Np. 

There are several combinatorial polynomial time algorithms for computing an integral 
maximum submodular flow, under the assumption that an oracle of computing the exchange 
capacity is available. They are designed by extending existing max-flow algorithms; see the 
survey m for further detail. 

We note one basic property for the case where the network is skew-symmetric. 

Lemma 2.2. Suppose that V is a signed set with t = s, and N is skew-symmetric. If p satisfies 
p{X_) < p{X) for all X CU, then the unique minimal minimum {s,t)-cut is a transversal. 

Proof. By the assumption for p, it suffices to show c{6lQ < c{5X). Suppose that some edge 
uv{^ dX) appears in 6X. In this case, it holds u G X^F X and u G X \ X. This means that 
{u, u} C X and n 0 X. Hence edge vu G 5X (of the same capacity) does not appear in dX. 
Consequently the cut-capacity does not increase. □ 
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2.2 Submodular extension 

Let V := {1, 2,..., n}. Let 3^ denote the set of ordered pairs of disjoint subsets of V. For a 
function h on 3^^, the polyhedron 'D{h) in is defined to be the set of all vectors z G R'^ 
satisfying 

z{Y) - z{Z) < h{Y, Z) ((y,Z)G3^). (2.2) 

If /i is a bisubmodular functiorQ then 'D{h) is known as a bisubmodular polyhedron. We study 
the case where T){h) is a projection of the base polyhedron of some submodular function. 
A representative of such polyhedra is the polyhedron of all flow-boundaries of a bidirected 
network; see [1]. 

Let = {I"*", 2+,..., n"*", 1“, 2“,..., n“} be the signed extension of V. For a function h 
on 3^^, a normal submodular extension of /i is a submodular function p on 2^"*" with p{%) = 0 
satisfying 


p{Y+yjZ-) = h{Y,Z) {iY,Z)€3^), (2.3) 

p(X) = p{X) < p{X) (Ag2^*). (2.4) 

Define a map (projection) (]) : R'^* —)■ R^ by 

(<^(^))(i) := (^eRV^±,iey). (2.5) 

Lemma 2.3. Let h be a function on 3^^ and p a normal submodular extension of h. Then it 
holds <f{B{p)) =V{h). 

Proof. Note that p{V^) = p{%) = 0. We first show that (j){B{p)) C 'D{h). Take an arbitrary 
X in B{p). For {Y, Z) G 3^, let A := y+ U Z”. Note that A = A and p(A) = p{fQ = p{X). 
Then we have 

x(A) < p{X), x{X) < p{X) = piX), x(y±) = p{V^) = 0. 

Hence we have (x(A) + x(A) - x{V^))/2 < {p{X) + p(X) - p{V^))/2 = p{X) = h{Y,Z). 
Also we have 

(x(A) + x(A) - x(y^))/2 = x(y+) + x{z-) - (x(y+ u z+) + x(y“ u Z"))/2 

= (x(y+) - x{Y-))/2 - (x(Z+) - x(Z“))/2 

= (<(-(x))(y)-(<^(x))(z). 

Hence (f{x) belongs to V{h). 

Next we show the converse. Take an arbitrary 2 : G D(/i). Define a vector x in R^"*" by 
x{i'^) := z{i) and x{i~) := —z{i) for i € V. Obviously 4>{x) = z. It suffices to show that x 
belongs to B(p). Since x{V^) = 0, we have x{V^) = 0 = /9(0) = p{V^). For A C V^, we 
have 

x{x) = x{x)= Y, E = <KA). 

i:i+eX 

Thus X G B{p), and hence D(/i) C (p{B{p)). □ 

^ A function h on 3^ is called bisubmodular if it satisfies 
h{X,Y) + hiX',Y') > h{XnX',YnY') + h{{XuX')\{YuY'),{YuY')\{XuX')) ((A, T), (A', A) G 3'^). 
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If /i : 3^^ —)• R has a normal submodular extension, then h is necessarily a bisubmodular 
function. Not all bisubmodular functions admit submodular extensions (Y. Iwamasa 2015). 

We consider a special bisubmodular function on 3-element set {1,2,3}, which plays a key 
role in Section]^ For 6 > 0, let A^, be the function on 3 ( 1 . 2 , 3 } (defined by 

( 2b if \Y\ > 2, 

A,iY, Z) :=\ b if |Y| = 1, |Z| < 1, (Y, Z) G 3^. (2.6) 

[ 0 otherwise (Y = 0 or \Z\ > 2), 

f 1 2 3 ! 

Lemma 2.4. The polyhedron T){Ai,) is the set of nonnegative veetors z G R)|_ ’ ’ ‘ satisfying 

z{l) + z{2) + z{3) < 2b, 
z{l) - z{2) - z(3) < 0, 

-z(l) + z(2) - zi3) < 0, 

-z{l) - zi2) + z{3) < 0. 


Proof. Observe that these inequalities appears in 


So it suffices to show that inequali¬ 


ties (2.2) are implied by the above inequalities. This is a routine verification. For example, 
z{l) + z{2) — z{3) < Z\;,({1, 2}, {3}) = 2b is obtained by adding z{l) + z{2) + z{3) < 2b 
and —2;(3) < 0. Also z{\) < Z\f,({l}) = 6 is implied by z{\) + z{2) 2;(3) < 2b and 
z{l) — z{2) — z{3) < 0, and z{l) — z{2) < Z\;,({1}, {2}) = 6 is implied by z{l) — z{2) — z{3) < 0 
and ^(3) < b. □ 


The polyhedron V{Aif) is a simplex of vertices (0,0,0), (6,6,0), (6,0,6) and (0,6,6). We 
will see in Section that ^{Af) represents the flow-conservation law and the node-capacity 
constraint on a node of degree 3. This bisubmodular function A^, has a normal submodular 
extension. The following example was found by Yuni Iwamasa via a computer calculation. 
Notice that each X C {I'*', 2"*“, 3'*', 1“, 2“, 3“} is one of the following six types: 

type 1: > 2 and \X~\ < 1. 

type 2: X~^ = and X~ = {!“, 2“, 3“} \ {Yj for some i G (1, 2,3}. 
type 3: A C {l-,2-,3-} or {l-,2-,3-} C A. 
type 4: |A+| = 2 and |A“| = 2. 

type 5: A+ = \X~\ < 2, and X~ / {!“, 2“, 3“} \ {*“} for some i G {1, 2, 3}. 

type 6: A = (I"'', 2+, 3''', 1“, 2“, 3“} \ {i“} for some i G {1, 2, 3}. 

Define Z\^ : 2'^^"'"’^''^’^^’^ ^ R by 


26 if A: type 1, 

Al{X) 0 if A: type 2 or 3, 

6 otherwise (A: type 4, 5, or 6). 

Lemma 2.5. AT is a normal submodular extension of Af,. 


(2.7) 


Proof. It suffices to con sider the case of 6 = 1; we denote Z\i and A\ by A and Z\*, respectively. 

First we show (2.3). For {Y, Z) G 3^^’^’^^, let A := Y+ U Z~. Then A+ = Y+ and 
X~ = Z~. If |Y| = lA+l > 2, then \Z\ = \X~\ < 1, and A is of type 1; hence A*{X) = 
2 = A{Y,Z). If |Y| = |A+| = 1 and \Z\ = |A-| < 1. Then A is of type 5, and hence 
A*{X) = 1 = A{Y,Z). If Y is empty, then A C {1“,2“,3“}, and A is type 3; hence 
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A*{X) = 0 = A{Y, Z). li\Z\ = \X-\>2 and Y then |y| = |X+| = 1, and X is of type 
2; hence A*{X) = 0 = A{Y, Z). 

Second we show (2.4). It suffices to show that A*{X) = A*{X) holds for any transversal 
X, and that A*(JQ < A*{X) holds for any X that is not transversal and not co-transversal. 
The former property follows from = 1 = Z\*({T'', 2''', S"'', 2“, 3“}), Z\*({1“}) = 

0 = Zi*({2+,3+,l-,2-,3-}), Z\*({l+,2+}) = 2 = 2+, 3+, 3-}), Zl*({l+,2-}) = 1 = 

Z\*({1^, 3''', 2“, 3“}), and Z\*({1“,2“}) = 0 = Z\*({3^, 1“, 2“, 3“}). The latter property 
follows from Zl*({2+,2-}) = 1 > 0 = Z\*(0), Z\*({1+, 2+, 2“}) = 2 > 1 = Z\*({1+}), and 
Z\*({2+,l-,2-}) = 1 > 0 = Z\*({1-}). 

Finally we show the submodularity of Z\*. Take X,Y C {!+, 2+, 3'’', 1“,2“,3“}. We can 
assume that X <^Y and T ^ X. 

Case 1: X is of type 6. In this case, X U X is the whole set, and is of type 3. Therefore it 
suffices to consider the case where X n X is of type 1 and X is not of type 1. Necessarily X 
is of type 4 or 6. Thus submodular inequality 1 -|- 1 > 2 + 0 holds. 

Case 2: X is of type 2 with I"*" G X. If X is also of type 2, then X U X(D {1~, 2“, 3“}) is 
of type 3, X n X(C {!“, 2“, 3“}) is also of type 3. If X is of type 3 with X C {!“, 2“, 3“}, 
then both X H X and X U X are of type 3. If X is of type 3 with X D {!“, 2“, 3“}, then X 
cannot have I"*", and thus both X n X and X U X are of type 3. In these cases, submodularity 
(0 + 0 > 0 -|- 0) holds. Thus we may assume that X is of type 1,4, or 5. Observe that neither 
X n X nor X U X is of type 1. We may assume that X contains 1"*“ and does not contain 1~. 
otherwise X n X or X U X is of type 3, and submodularity holds. Necessarily X is of type 1. 
Then X 0 X is of type 5, and X U X is of type 4 or 6; submodularity (0 + 2 > 1 + 1) holds. 

Case 3: X is of type 3 and X is not of type 2. If X is also of type 3, then both X n X 
and X U X is of type 3; submodularity holds. Since one of X n X and X U X is of type 3, 
it suffices to consider the case where X n X or X U X is of type 1. We show that X is also 
type 1. If X n X is of type 1, then X = {1+, 2+, 3''', 1“, 2“, 3“} \ {i^}, and necessarily X is 
of type 1. If X U X is of type 1, then X = {i~}, and X consists of at least two elements in 
{l+,2+,3+} (type 1). 

Case 4: X is of type 4 or 5, and X is type 1, 4, or 5. Suppose that X is of type 4. Then 
XUX is not type 1. We may consider the case where XnX is of type 1 and X is not of type 1. 
Then X contains X"*". Thus X is of type 4, and XU X necessarily contains {1“, 2“, 3“} (type 
3); submodularity (1 + 1 > 2 -|- 0) holds. Suppose that X is of type 5 with 1+ G X. Then 
X n X is not type 1. We may consider the case where X U X is of type 1. If X does not 
contain 1+, then X n X is of type 3; submodularity (1 + 1 > 2 -|- 0) holds. Thus |X+| > 2 
and |X“| = 0 or 1; X is of type 1. The intersection XnX is of type 5; thus 1 -|- 2 > 2 + 1 
holds. □ 


3 Node-capacitated multiflow 

In this section, we introduce a combinatorial duality theory, developed by m. for a class 
of node-capacitated multiflow problems. We consider the following multiflow problem. Now 
assume that network N also has a nonnegative edge-cost a : E ^ R+; so the network is a 
5-tuple (X, E, S, c, a). For a multiflow / = {V, A), the total flow-values on node i and edge e 
are denoted by /(i) := ieV{P) /(®) •= EpeP: e&E{P) KP)^ respectively. The 

cost a{f) is dehned by 

a(/) := 

e£E 

Next we define the value of a multiflow. A tree-embedding E = (F, {^sjses) is a pair of a tree 
r and a family {( 7 s}sgs of vertices of E indexed by terminal set S. The £-value V£{f) of a 
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multiflow / = (P, A) is defined by 

^sif) ■= ^(^^p’9tp)A(P), 

P€P 

where sp,tp denote the ends of an 5-path P, and d = dp denotes the shortest path metric 
of r with respect to unit edge-length. We are now ready to define our multiflow problem. 
An instance of the problem is a pair of a network N = {V, E, S, c, a) and a tree-embedding 
S = {r, {ps}se5)) and the task is to find a multiflow / that maximizes V£{f) — a{f). 

This somewhat artificial formulation turns out to be useful, and actually generalizes the 
original problem. Indeed, take T as a star with |5| leaves Vs (s G 5), let £ := (r,{us}se 5 ), 
and let a(e) := 0 for each edge e. Then V£{f) — a{f) is twice the total flow-value of /. 

In Section |3.1[ we present a combinatorial duality theorem and an optimality criterion. 
We introduce a nondegeneracy concept of the problem, and give an algorithm to find a half¬ 
integral optimal multiflow from a dual optimum under the nondegeneracy assumption. We 
also explain how to reduce the original problem to a nondegenerate problem. In Section [3.2t 
we see that our dual objective can be viewed as an L-convex function on a certain graph 
structure, and present the steepest descent algorithm (SDA) to minimize L-convex functions 
and its iteration bound. 


3.1 Duality 

Let a pair of A" = {V, E, S, c, a) and £ = {r,{qs}ses) be an instance of the problem. The 
vertex set of P is also denoted by P (instead of V{r)). Let P* denote the edge-subdivision 
of P, where P C T*, the edge-length of P* is defined as 1/2 uniformly, and the shortest path 
metric dp* is also denoted by d. 

A pair (p, r) of a tree-valued function p : V ^ P* and a nonnegative half-integer-valued 
function r : 1/ —)■ Z+/2 is called a potential if it satisfies the following conditions: 

(pi) For each node i, r{i) is an integer if and only if p{i) belongs to P. 

(p2) For each edge ij, it holds d{p{i),p{j)) — r{i) — r{j) < a{ij). 

(p3) For each terminal s, it holds (p(s),r(s)) = (gs,0). 


Then the following min-max formula and optimality criterion hold: 


Theorem 3.1 ([IS]). Suppose that a is even-integer-valued. The maximum of vg{f) — a{f) 
over all multiflows f is equal to the minimum o /potentials {p,r). 

Lemma 3.2 ([IS])- Suppose that a is even-integer-valued. A multiflow f = ('P,A) and a 
potential {p, r) are both optimal if and only if they satisfy the following conditions: 


(ol) For each path P in V with X{P) > 0, it holds d{qsp, qtp) = E d{p{i),p{j))- 

ij&E{P) 

(o2) For each edge ij with f{ij) > 0, it holds d{p(i),p(j)) — r{i) — r(j) = a{ij). 

(o3) For each nonterminal node i with r{i) > 0, it holds f{i) = c{i). 


We will use the if part (and the weak duality in Theorem 3.1) only, which is proved for 
completeness. 
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Proof. (If part). For any multiflow / = {V,X) and any potential {p,r), the difference 
Eiey\s2ciri - (veif) - a{f)) is equal to 

^ 2(ci - fi)ri + ^ fijiaij - d{pi,pj) + r* + rj) 
i£V\S ijeE 

+ Ew E d{pi,Pj) - d{qsp,qtp)\ >0, (3.1) 

PSP \ij&E{P) j 

where we use 


X] fiAPi^Pj) = Y. d{Pi,Pj), 

ij&E Pev ij&E{P) 

Y = Y + 

i£V\S ij&E 


Thus, if / and {p,r) satisfy conditions (ol), (o2), and (o3), then the equality holds in (3.1), 
and both / and {p, r) are optimal. □ 


Nondegenerate case. An instance {N,£) is said to be nondegenerate if the edge-costs a 
are positive even-valued and the degree of each node in P is at most 3. Suppose that {N,S) 
is nondegenerate. We further assume, for notational simplicity, that tree P has no vertex of 
degree one (by attaching paths of infinite length). Let r 2 and Ps denote the sets of vertices of 
P with degree 2 and 3, respectively. For a vertex v in P, the neighbors of u in F are denoted 
by u_s.i,u_s .2 if V G r 2 and u_ 5 .i, u_>. 2 , u-s-s if u G Ps- Let P* denote the edge-subdivision of P. 
For vertex v G P*, the neighbors of v in P* are denoted by v^*2 if v G r2 or v G P* \ P, 
and v^*i,v^* 2 ,v ^*3 ifvGPs- Let denote the connected component of F* — u containing 

We are going to characterize the flow support of an optimal multiflow. Let {p, r) be a 
potential. Motivated by (o2), define the edge subset Ep^r by 

Ep,r ■= {ij G E I d{p{i),p{j)) - r{i) - r{j) = aij}. 

For a nonterminal node i, let 5p,fc(0 denote the set of edges ij G Ep^r with p{j) G P^ fc' 
each edge cost a{ij) is positive, it holds p{i) / p{j) for ij G Ep^r- Thus 5p^k{P) for A: = 1,2,3 
(or k = 1,2) partition the set (5{z} of all edges in Ep^r incident to i. 

A {p, r)-admissible support is a function ( : Ep^r —^ IL+ satisfying the following conditions: 

(al) For each nonterminal node i with pii) 0 F 3 , it holds C(<5p,i(*)) = C{dpfi{i)) < c(i). 

(a2) For each nonterminal node i with p{i) G F 3 , it holds 

C,{dp,i{i)) + C(dp, 2 (i)) + C(dp,3(i}) < 2 c(i), 

C(<Jp,i(*)) - C(dp,2(i)) - C(dp,3(l)) < 0 , 

-C('^P,i(*)) + C(V2(*)) - CA,3(^)) < 0, 

-C(Vi(*)) - C(dp, 2 (i)) + C(dp,3(^)) < 0 . 

(a3) For each nonterminal node i with r(i) > 0, it holds C(<^{*}) = 2c(i). 

(a4) For each edge e, C(e) is a half-integer, and for each nonterminal node i, C(d{i}) is an 
integer. 
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Recall the notational convention C{^p,k{i)) ■= ^ (p, r)-adniissible 

support is viewed as an edge-weight whose degree vector belongs to a bisubmodular polyhedron 
described by /ic(i)- 

It is not difficult to see from Lemma |3.2| that for any half-integral optimal multiflow /, 
the flow-support C of /, defined by C(e) := /(e), is a (p, r)-admissible support. The converse 
also holds. 


Lemma 3.3 ([E]). Let {p,r) be a potential. If a {p,r)-admissible support C exists and is 
given, then {p,r) is optimal and a half-integral optimal multiflow is obtained in 0{nm) time. 


Thus our problem is to find a potential {p, r) so that a (p, r)-admissible support exists. 
An algorithm for Lemma 3.3 is the following. 


Algorithm 1: Construction of an optimal multiflow from a (p, r)-admissible support. 


Input: A potential (p, r) and a (p, r)-admissible support (. 
Output: A half-integral optimal multiflow / = {V, A). 


Step 0: V = %. 


Step 1 : Choose a terminal s and an edge sj with C('Sj) > 0- If such a terminal does not 

exist, then / = {'P,\) is a half-integral optimal multiflow; stop. Otherwise let jo ■<— s, 

Ji ^ j, C{sj), I ^ 1, and go to step 2. 

Step 2: If ji is a terminal, then add path P = ... ,ji) to P with flow-value A(P) := p, 

let C(e) •(— C(e) — p for each edge e in P, and go to step 1. Otherwise go to step 3. 

Step 3: If p{ji) 0 P 3 and ji-iji G dp^kUi) for k G {1,2}, then choose an edge jiji+i from 

Sp,k'Ui) with k' k and C(j)//+i) > 0, and let p ^ minjp, C(/ij}+i)}. 

If pUi) ^ -I 3 and ji_iji G 5p,fc(//) for k G {1,2,3}, then choose an edge ji+iji from 
^p,k'{jl) with k' k, C{jiji+i) > 0 , and C{^p,k{jl)) + C{Sp,k’{jl)) - C{Sp,k"{jl)) > 0 for 
fc" G {l,2,3}\{fe,A:'}. Let 


p min 




C{^p,kUl)) + C{^p,k'{jl)) - C{^p,k"ijl)) I 


Let I I + 1 and go to step 2. 


This algorithm is essentially the proof of |15[ Lemma 4.5]. Let us sketch the correctness of the 
algorithm; we show that the resulting multiflow / satisfies the conditions (ol),(o2), and (o3) 
with (p, r). The condition (o2) follows from /(e) = 0 for e G E\Ep^r- In step 3, 


3.2 


in Lemma 

we can always choose a required edge by (al) and (a 2 ). Also C still satishes the conditions 
(al), (a2), and (a4), thanks to the way of the update. Each produced path (/o,Ji, J 2 , ■ • • ,jm) 
satisfies 


d{p{ji-i),p{ji)) + d{p{ji),p{ji+i)) = d{p{ji_i),p{ji+i)) (1 < / < m - 1) 

since p(j/-i) G PUi-ei) ^ ^p{ji),k' ^ ^ d{p{ji-i),p{ji)) is positive 

(since a is positive). Since T is a tree, we have d{p{jo),p{jm)) = YdLid{p{ji-i),p{ji))\ 
see e.g., iia Lemma 3.9]. Thus each produced path satishes (ol). By the same argument, 
every edge e with (j{e) > 0 extends to an S'-path consisting of edges e' with C,{e') > 0 
satisfying (ol). This means that if no terminal s is chosen in step 1, then (} = 0. By (a3), 
the resulting multihow / satishes (o3). Notice that p is a half-integer by (a4). Hence / is 
half-integral and optimal. Once an 5-path P is obtained, Q becomes zero on some edge, or 
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C(Vfc(^))+C('5p,fc'(i))-C(<5p, k"{i)) becomes zero on some node i; they remain zero in subsequent 
iterations. Thus the algorithm terminates after 0{m) paths are obtained, where each path is 
found in 0{n) time by keeping {e G Sp^k{i) \ C(e) >0} {i G V, k = 1, 2, 3) as lists. 

We estimate the range in which an optimal potential exists. Let Fq denote the minimal 
subtree in F containing {gsjsgs, and let fi(T'o) denote the diameter of Tq, i.e., (i(T'o) := 
max^^^g^o d{u,v). 

Lemma 3.4. There is an optimal potential {p,r) with p{i) G Fq and r{i) < d(T'o) for i gV . 

Proof. Let (p, r) be a potential. Suppose that there is a nonterminal node i* with pi* 0 Fq. 
Then there is s G S' such that the path between qg and pt* is outside of Tq) Take such s, i* with 
d{qs,Pi*) maximum. We can assume that contains qg. Let X be the set of nodes j with 

Pj = Pi*- Suppose that pi* G F* \ F. Then rj > 1/2 for all j G X. For each j G X, replace 
(PjXj) by {pj^*i,rj — 1/2). For an edge ij with i G X and j ^ X, both d{pi,pj) and r* + rj 
decrease by 1/2, and thus (p2) remains to hold. For other edge ij, quantity d{pi,pj) —ri — rj is 
nonincreasing or remains nonpositive (if i,j G X). The feasibility (p2) still holds (since a{ij) 
is nonnegative). Thus the resulting {p,r) is a potential, and the objective value decreases. 
Suppose that pi* G F. For each j G X, replace {pj,rj) by {pj^.,^,rj). For each edge ij, 
distance d{pi,pj) does not increase. Thus the feasibility (p2) holds, and the objective value 
does not change. By repeating this procedure, we can make (p, r) so that pi G Fq for i G V, 
without increasing the objective value. Suppose that r* > d(/o) for some i; necessarily > 1. 
For each edge ij connecting i, d{pi,pj) — ri — rj — aij < —1 (since d{pi,pj) < d{FQ) and 
d{pi,pj) — ri — rj is an integer). Thus we can replace r, by — 1 to decrease the objective 
value. Repeating this procedure, (p, r) satisfies r* < d{Fo), as required. □ 

Reduction to a nondegenerate instance. Here we explain how to reduce our original 
problem to a nondegenerate problem. An instance of the original problem is viewed as a pair 
of network N = {V, E, S,c,a) and a tree-embedding S = {F,{vg}s^s) such that a(e) = 0 
for all edges e and T is a star with center vq and leaves Vg (s G S). We are going to 
construct a nondegenerate instance. Define edge cost a by a(e) := 2 for each edge e G E. 
Let N := (y,E,S,c,d). Next we define a tree-embedding S = {F,{qs}g^s)- Let E be any 
(finite) trivalent tree with l^l leaves Ug {s G S) and diameter D = 0(log IS*!). For each s G S, 
consider an infinite path Pg having a vertex of degree one. Identify Ug and u^, i.e., glue Pg 
and S at Ug. The resulting infinite tree is denoted by F. Define qg as the vertex in Pg having 
distance {2\E\ + 1)D from Ug{= u^). See Figure [l| for the construction of F. 

Now we obtain a nondegenerate instance {N,S). Let (p, r) and / = {V,X) be an optimal 
potential and an optimal multiflow, respectively, for this perturbed instance {N,£). We show 
that / is a maximum multiflow, i.e., optimal for the original instance {N,£). We are going 
to construct an optimal potential (p, r) for {N,£) from {p,r). Let Bi denote the subgraph of 
F induced by vertices q with d[pi,q) < fi. Namely Bi is (the intersection of F and) the ball 
with center pi and radius fj. Then it holds 

d{Bi,Bj) := min d{u,v) = mm{0, d{pi, pj) — hi — hj}. (3.2) 

u£Bi,vGBj 

For k = 0,1,2,..., 2\E\ and s G S, let Fg^k denote the subgraph of F consisting of edges uv 
in Pg C F such that kD < d{ug,u) = d{us,v) — 1 < (fe -|- 1)D. Let Fk be the union of Fg^k 
over s G S. We say that an edge ij in E hits Fk if Bi n Rj = 0 and the path between Bi and 
Bj meets an edge of Fk. For an edge ij, it holds d{Bi,Bj) < Oij < 2. Therefore the path 
between Bi and Bj consists of at most two edges. Since 2|R| -|- 1 subgraphs Fq, Fi, ..., F 2 \e\ 
are edge-disjoint, there is an index k such that every edge in E does not hit Fk. Fix such an 
index k. 
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Figure 1: Construction of F 


For each s £ S, choose the edge of Fg^k nearest to qg (furthest from Ug). Delete all eg 
from F. There are |5| + 1 connected components Co,Cg {s G S), where Cq is the connected 
component containing S, and Cg is the connected component containing Qg. For each i £ V, 
define {pi,ri) £ F* MZ* by 

{ (uo, 1) if Bi contains two e*, eg', 

(us,l/2) if Bi contains exactly one eg, (3-3) 

{vg, 0) if Bi is contained in for s G 5 U {0}, 

where Vg denotes the vertex in F* \ F obtained by subdividing edge VQVg in F. We show that 
{p,r) is a potential for {N,S) and satisfies (ol), (o2), and (o3) with /. We first show the 
feasibility (p2) dr{pi,Pj) < + rj for ij £ E. Since pi £ F* \ F implies r* = 1/2, we may 
consider the three cases: (i) pi = Vg, pj = vq, (ii) pi = Vg, pj = Vg' for s / s' , and (iii) pi = Vg, 
Pj = Vg' for s ^ s'. For (i), Bj necessarily contains two eg, eg'-, otherwise ij hits F^ at e*. This 
implies Vj = 1. (ii) and (iii) cannot occur since, otherwise, ij hits F^ at eg. Thus (p2) holds, 
and hence {p, r) is a potential; (pi) and (p3) are clearly satisfied. 

Next we show (ol),(o2), and (o3) for / and {p,r). To show (o2), take an edge ij with 
f{ij) > 0. By (o2) for / and {p,f) in {N,£), it holds dp{pi,pj) — Vi — fj = 2. Thus the balls 
Bi and Bj are disjoint and have distance 2. Suppose that Bi has two edges eg, eg'. Then Bi 
must meet Fg^^ for every s £ S. Necessarily Bj cannot be contained by Cq. Also Bj cannot 
have et for any t £ S\{s, s'}; otherwise ij hits F^. Thus Bj is contained in Ct for t £ S, and 
{Pi,ri) = (uo, 1) and {pj,rj) = {vt,0) hold. If Bi has (only one) eg and Bj has (only one) e^/, 
then s and s' must be different, and necessarily {pi,ri) = {vg,l/2) and {pj,rj) = {vs',1/2). 
If Bi has only eg and Bj does not have any of et, then necessarily Bj is contained in Cg or 
Co; hence {pi,ri) = {vg,l/2) and {pj,rj) = (^ 5 , 0 ) or (uo,0). If both Bi and Bj do not have 
any of eg, then both Bi and Bj are contained in Cg for some s G S' U {0}, and pi = pj and 
Ti = Tj = 0. In all the cases, it holds dr{pi,Pj) = r* + Vj, implying (o2). 

Consider the condition (ol). Take a path F = (s = jo,ji,... ,jm = t) with A(F) > 0. 
There is an index I such that Bi contains e^; otherwise F^ is hit by some edge. Moreover 
such an index I is unique. Otherwise, the balls Bj^ and Bj^ with I < I' contain eg. Then 
diPji,Pji,) — fji — fji, < 0. However, by (ol) and (o2) for / and {p,r), we have d{pj^,pj^,) = 

Oi + Oi+i + 2, and d{pj^,pj^,) - fj^ - fj^, > 2 > 0; a contradiction. 
Similarly there is a unique index V such that Bf contains e*. If / = /', then {pjQ,Pji,.. ■ ,Pj„^) 
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Figure 2: Graph F* Kl Z* 


must be {vs,Vs, ■ ■ ■ ,Vs,vo,vt,... ,vt). If I / I', say I < I', then {Pjo,Pji, ■ ■ ■ ,Pjrr,) must be 
(vs,. . ■ ,Vs,Vs,vo, .. . ,vo,vt,vt, ... ,vt) or {vs,... ,Vs,Vs,vt,vt,... ,vt). Thus we obtain (ol). 
Since f(i) = 0 implies r{i) = 0, we obtain (o3). Hence {p, r) is an optimal potential, and / is 
a maximum multiflow. 

3.2 Discrete convexity and steepest descent algorithm (SDA) 

Here we briefly introduce a class of discrete convex functions (L-convex functions) on a certain 
graph structure and the steepest descent algorithm to minimize them. We then explain that 
our problem falls into the minimization of an L-convex function. The full discussion will be 
given in |18j . 

First we equip the space of all potentials with a graph structure. Let Z*(;= Z/2) denote 
the set of half-integers. Let F* Kl Z* denote the set of pairs {p, r) G F* x Z* such that p G F 
if and only if r G Z. Two points {p,r) and {p',r') are adjacent if and only if p and p' are 
adjacent in F* and |r — r'| = 1/2. Fix an arbitrary vertex po of F. Let B (resp. W) denote 
the subset of F* Kl Z* consisting of pairs {p, r) with p G F, r G Z, and d{p,po) + r even (resp. 
odd). Orient each edge of F* Kl Z* by {p,r) G- {p\r') if (p,r) G H or {p',r') G W. Namely 
B is the set of sinks and W is the set of sources. This orientation is acyclic, and induces a 
partial order ^ on F* Kl Z*. See Figure where nodes in B and W are colored black and 
white, respectively. 

Next we define midpoint operations on F* MZ*. Let F** denote the edge-subdivision of 
F* with edge-length 1/4, let Z**(;= Z/4) denote the set of quarter-integers, and let F** Kl Z** 
denote the set of pairs {p,r) G F** x Z** such that p G F* if and only if r G Z*. For two 
points X = (p, r),x' = (p', r') in F* Kl Z*, there exists a unique midpoint y = {q, t) G F** Kl Z** 
such that d{p, q) + d{q,p') = d{p,p'), d{p, q) = d{q,p'), and t = {r + r')/2. This y is denoted by 
(x + x')f2; accordingly q is denoted by (p -|- p')j2. For z = (g, t) G F** Kl Z**, there uniquely 
exists a pair (x, y) of vertices in F* lEI Z* with the property that z = {x + y)/2 and x <y. We 
denote x and y by [z\ and [ z ], respectively. 
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We are ready to define L-convex functions. For a natural number n, consider the product 
(r* lEI Z*)”; a point x in {F* Kl Z*)^ is represented by a pair {p, r) of p £ and r G (Z*)”. 

A function g : (F* Kl Z*)"■—)■ R U { 00 } is L-convex if it satisfies the following analogue of the 
discrete midpoint convexity [221 Section 7.2]: 

g{x)+g{y)>g{[{x + y)/2\) + g{\{x + y)/2\) {x,y £ {F* MZ*)^), (3.4) 

where {[{x + y)/2\)i := [{xi + yi)/2]\ and {\{x + y)/2~\)i := \{xi + yi)/2~\ for i = 1, 2 ,..., n. 

For X £ {F* Kl Z*)”, let iFx (resp. Ix) denote the set of points y with Xi ■< yi (resp. 
Xi F yi) for i = 1, 2,..., n. The set Fx U Xx is called the neighborhood of x. The steepest 
descent algorithm is given as follows: 

Algorithm 2: Steepest descent algorithm (SDA) 

Input: An L-convex function g : [F* Kl Z*)"' —> RU { 00 }, and a point x^ in domgf. 

Output: A minimizer of g. 


Step 0: Let i ■<— 0. 

Step 1 : Find a minimizer y of g over the neighborhood Fxi UX^i of x*. 
Step 2: If g{x^) = g{y), then output x* and stop; x* is a minimizer. 


Step 3 : Otherwise, let x 




y, z •(— i + 1 , and go to step 1 . 


The fact that the output is a minimizer easily follows from (3.4); see nil Theorem 2.5]. 

We discuss the number of iterations of this algorithm. For x,y £ (F* lEI Z*)^, an loo-path 
between x and y is a sequence P = {x = x°,x^,... ,x^ = y) such that for each k and i, 
the i-th components x^ and x^^^ belong to a common 4-cycle of F* Kl Z*, where m is called 
the length of P. The loo-distance between x and y, denoted by Doo{x,y), is defined as the 
minimum length of an / 00 -path between x and y. For distinct {p,r), {q, s) G X* Kl Z* in a 
4-cycle, it holds d{p,q) -|- |r — s| = 1. From this we observe that 


Doo{x, y) = max d{pi, qi) -h \ri - Si 

l<i<n 


ix = {p,r),y = {q,s) £ (X* K Z*)"). 


(3.5) 


Notice that a sequence (x = x'^jX^jX^,... ,x”^) generated by SDA is an ^oo-path. Let opt{g) 
denote the set of all minimizers of g. The length m, which is the number of the iterations, is 
at least Xoo(xo, opt( 5 )) = minygopt(g) Dooixo,y)- This lower bound is almost tight. 

Theorem 3.5 ([IS]). Let (x = x°, x^,..., x™) be a sequence of points generated by SDA 
applied to an L-convex function g and an initial point x. Then m < Doo{x,opf{g)) -\- 2. If 
g{x) = minygjy^ g{y) or g{x) = minygi^ g{y), then m = Xoo(x, opt(c/)). 

A similar bound for original L-convex functions in DCA was established in [23] . For a 
similar but different class of L-convex functions (called alternating L-convex functions), the 


same bound was proved by m- By using this result, we give a shorter proof of Theorem |3.5| 
in Appendix. 

We now return to our problem. Let a pair of A = {V, E, S, c, a) and £ = (X, {(?s}ses) be 
an instance. Let V = {l,2,...,n}. Then the set of potentials (p,r) £ (F*)^ x (Z+/2)^ is 
naturally regarded as a subset of (X* lEI Z*)”". Define a function gj\f^£ : (F* Kl Z*)*^ — )> RU {00} 
by 

gN,s(x) ■■= ^ 2c(i)r(i) (3.6) 

iev\s 

if X = (p,r) is a potential and gN,£(x) '■= 00 otherwise. 

Propositiou 3.6 ([IS]). Suppose that a is even-valued. Then gN,£ is L-convex. 

In Appendix, we give a shorter proof of this proposition for our case where a is positive. 
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4 Algorithm 


In this section, we prove the main result (Theorem 1.1) by presenting an algorithm {dual 
descent algorithm) to solve a nondegenerate instance of our multiflow problem. We first 
show that the optimality check of a potential {p,r), or finding a (p, r)-admissible support, 
is reduced to the submodular flow feasibility problem on a certain skew-symmetric network, 
called the double covering network. This extends the earlier result on the minimum-cost edge- 
capacitated multiflow problem by Karzanov [191 [20] , in which the optimality is checked by the 
classical circulation problem. Partial adaptations of this idea to the node-capacitated setting 
have been given in [am; but the full adaptation using submodular flow is new. Checking the 
feasibility of a submodular flow is reduced to the maximum submodular flow problem. We 
prove that the minimal minimum cut naturally gives a steepest direction at each potential. 
Then we obtain a simple descent algorithm mentioned in Introduction. 


4.1 Double covering network with submodular constraints 

Let a pair oi N = {V,E,S,c,a) and S = {F, {qs}s£s) be a nondegenerate instance. We 
can assume that there is no edge such that both of its ends are terminals. Let (p, r) a 
potential for {N,S). A nonterminal node i is said to be flat if p{i) 0 T 3 , and singular if 
p{i) G r^. A nonterminal node i is said to be zero if r{i) = 0 and positive if r{i) > 0. 
We are going to construct a skew-symmetric network associated with (p, r). We first 
define a signed node set Ui and edge set Ai indexed by each node i in N, together with 
lower edge-capacity c and upper edge-capacity c. For each terminal s, let Ug := {s''',s“}, 
and let Ag := with c(s^s“) := 0 and c(s'''s“) := 00 . For each flat node i, let Ui := 

{i^, Z 2 }, and let Ai := {i'^ with = c(z^z^) := c{i). If i is positive, then 

cfli'^F^) = cflii^ii) '■= c{i). Otherwise := 0. For each positive singular node 

i, let Ui := {iq fl'l,i 2 utU q fli ^2 A consist of i^i^, = 1)2,3) 

with = c{i^io) := 2c(i), c(i^io) = c(zoZ^) := 0, and cjz^z+)^= ^i^i'^) := c{{) 

{k = 1,2,3). For each zero singular node i, let Ui := {z^, Zg", z),", and let Ai := 0. 

Define a submodular function A* on Ui by A*{X) := where X' is obtained from 

X by replacing each z^ (resp. z^) with k'^ (resp. k~). For an edge e = ij £ Ep^r with 
ij e <lp,fc(*) and ij G Sp^k'iJ), add edges e+ = z^j^ and e“ = j^,i^ with c(e+) = c(e“) := 0 
and c(e''') = c(e“) := 00 . In the case where j is a terminal s, let = z^s^ and e~ = s~i^. 
Now the double covering network Fp^r = (U, A, c, c) is defined by the disjoint union of all these 
edges and nodes. See Figure [^ where z, j, z', and j' are positive singular, zero singular, zero 
flat, and positive flat nodes, respectively. 

A node v in Fp^r is said to be unusual if v belongs to U for a zero singular node z. Other 
nodes v are said to be usual. 

A circulation p in Fp^r is a function on edge set A such that 

c(e) < (p(e) < c(e) (e G A), 

(V<p)(z;) = 0 (usual node v in Fp^r), 

{Xip)\ui G B{A*) (zero singular node z in N). 


For a circulation p, define : Ep^r —^ H,+ by 

(4. 

Lemma 4.1. Suppose that an integral eirculation p in Fp^r exists. Then Cip zs a {p,r)- 
admissible support. 


Cv{e) := 


ip{e+) + ip{e ) 


(e G Epr). 
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Figure 3: Double covering network 


Proof. The half-integrality of (in (a4)) is clear. Take a nonterminal node i. Let fk ■= 
CipiSpM^)) = Eee<5j,.fe(i)(</^(e’^)+¥’(e"))/2. Suppose that i is flat. Then Ci = Eee5p.i(i)(<^(e’^) + 
</^(e“))/2 = {ipiifif) + (p{ith))/2 = Eee5p.2w(‘i^(®’^) += ^ 2 - Since (filitif) < c{i) 
and tp{i 2 if) < c{i). We obtain (al). In addition, if i is positive, then ip{ifif) = = 

c(i) must hold, and we obtain (a3). Since Ci and C 2 are half-integers with Ci = C 2 ) we have (a4). 

Suppose that i is positive singular. Since Ck = Yle&Spf,(i) {ip{e+) + ip{e ))/2 = {ip{i+i+) + 
V?(*o *fc))/2 < c{i), we obtain Cfc < c(i). Moreover it holds = C1 + C2 + C3 = <y^(*d^o ) = 

2c(i), implying (a3) and (a4). If Ci > C2 + Cs; then Ci + C2 + Cs < 2Ci < 2c(i); this contradicts 
Cl + C2 + Cs = 2c(i). Therefore Ci, C2, and Cs satisfy (a2). 

Suppose that i is zero singular. Notice that (V(/7)(zd) = Eee(5 k{i) ~ 

-Eee<5p,fem‘^(®“)- (Ci,C2,C3) = (pi{Vp)\ui); see ([^ for Since (V<y9)|r/, e B{A*), 


because of Lemmas 2.3 and 2.5, the vector (Ci)C 2 ,C 3 ) satisfies (a 2 ). Since 0 = {V(p){Ui) = 


Efc=i, 2 , 3 (^¥^)(^fc) + Efc=i, 2 . 3 (^<^)(*fc)> we have Efc=i, 2 , 3 (W)(^) = “ Efc=i, 2 , 3 (^<^)(*fc ) and 
obtain (a4) by 

C^W}) = Ci + C2 + C3 = ^ (v</^)(i+) - (V^)(z^) = (V(^)(i+)GZ. 

^ fc=l,2,3 fc=l,2,3 


□ 


4.2 Dual descent algorithm: implementing SDA by submodular flow 

To check the existence of a circulation in 'Dp^r, we construct an instance of the maximum 
submodular flow problem. Add a super source and super sink a~. For each edge e = v~^u~ 
in Vp^r having nonzero lower capacity c(e) > 0 , replace v~^u~ by two edges v~^a~ and o+tt" 
with (upper) capacity c(e) (and lower capacity 0). Those edges are z^Zq positive singular 
nodes z and if if ,12 if for positive flat nodes z. The resulting (skew-symmetric) network is 
denoted by Vp^r, where modified edge sets are denoted by Ai (z G V) and the (upper) edge- 
capacity is denoted by c. Consider the maximum (a"'', a“)-submodular flow problem on Vp^r, 
where submodular function p on [/ is given as 

p{X):= Y. ^*{XAUi) {XCU). (4.2) 

2 : zero singular 
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Lemma 4.2. If {a^} is a minimum {a^,a )-cut inDp^r, then a circulation (p in T>p^r exists, 
and is obtained from any maximum flow p' by the following procedure: 

(A) For each edge e = v~^u~ in Vp^r having nonzero lower capacity, let p{e) := p'{a'^u~), 
and for other edge e in Vp^r, let p{e) := p'{e). 

Indeed, in this case, any max-flow saturates all edges leaving a'^ and all edges entering a“, 
and consequently the resulting p satisfies p{e) = c(e) = c(e) for all edges e having nonzero 
lower-capacity, and is a circulation of 

Next we show that the minimal minimum cut gives rise to a steepest descent direction of 
gN,£ at An (a'*', a“)-cut X is said to be normal if it satisfies the following conditions: 

(cO) X does not meet {s"*", for any terminal s € S. 

(cl) A is a transversal. 

(c2) For each positive node i, X D is equal to 0, or {f^}, and X n Ufl is equal to 0, 
U-,ovU-\{^t}. 

(c3) For each zero node i, X nUi is equal to 0, or } U U~ \ {\}- 

Then the following holds; the proof is given in the end of this subsection. 

Lemma 4.3. A normal minimum (a"*", a~)-cut exists, and is obtained from the unique minimal 
minimum cut X by applying the following procedure: 

(B) For each singular node i, if \X n Ui^\ = 2, then replace X by X U . 

Let Ux be the subset of U consisting of i'lfflfi for i ^ V \ S and k G {1,2,3} such that 
{p{i),r{i)) G W, or p{i) G F* \ F and (j){i)^*k,r{i) — 1/2) G W. Let Ujr be the node subset 
defined by replacing W with B in the definition of Ux- A normal cut X is said to be IF-normal 
if A n Lx = 01 and is said to be of I-normal if A n Ujr = 0. For an A-normal or X-normal 
cut A, define {p,r)^ by 


ip,r)^{i) 


(p(i)^.fc,r(i)-h 1/2) 

{p{i)^*kU{i) - 1/2) 
(p(f)^fc,r(z)) 
{p{i),r{i) + 1) 
{p{i),r{i) - 1) 

. (p(f),r(f)) 


ii Xf^Ui = {i+}, 

if AnL, = L-\{i^}, 

if AnL, = {i+}ULr\{i-}, 

if An Li = L+, 

if An Li = L-, 

if A n Li = 0 


(4.3) 


for each nonterminal node i G V \ S and {p,r)^{s) := {qs,0) for each terminal s G S. See 
Figure]^ for the correspondence between {p,r)^{i) and L n A, where nodes in Li n A are 
colored black. 


Proposition 4.4. Let X be a normal minimum (a+,a )-cut in Fp^r, o-nd let Xjr := A n Ujr 
and Ax := AnLx- Then {p, r)^^ is a minimizer of gN,£ over IFp^r, and {p, r)^^ is a minimizer 

ofgN,£ overXp^r- 

The proof is given in the end. We are now ready to describe our algorithm. 

Algorithm 3: Dual descent algorithm (for a nondegenerate instance) 

Input: A nondegenerate instance N = {V,E,S,c,a), E = {F, {qs}s£s)- 
Output: A half-integral optimal multiflow /. 
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Figure 4: Neighborhood and normal cut 


Step 0: Choose any vertex v in Fq, and let {p,r) be a potential defined by {p{i),r{i)) := 
{v, d{ro)) for i G 14 \ S' and (p(s), r(s)) := {qs, 0) for s £ S. 

Step 1 : Construct network Obtain an integral maximum {a~^,a )-flow p' and the 

minimal minimum (a^,a“)-cut X in T>p^r by a submodular flow algorithm. Make X 
normal by Procedure (B). 


Step 2: If X = {o"*"}, then {p,r) is optimal, obtain a feasible flow p in Fp^r from p' by 
Procedure (A), obtain a (p, r)-admissible support Ctf by (4.1), and obtain a half-integral 
optimal multiflow / from by Algorithm 1; stop. 

Step 3: Choose Y G {Xj.Xx) with gN,e{{p,rV) = ^'^H9N,£{{p,r)^^), 9N,£{{p,r)^^)}- 
Let {p,r) <— {p,r)^, and go to step 1. 


Theorem 4.5. The dual descent algorithm runs in 0{d{rQ)MSF{n,m,l)) time. 


Proof. By Proposition 4.4, (p, r)^ is a steepest direction of gN^g at (p, r). Hence the dual 


descent algorithm is viewed as the steepest descent algorithm applied to gN,£ that is L-convex 
(Proposition 3.6). By Lemma 3.4 and Theorem 3.5 the number of iteration is bounded 


by d(Po)- Our submodular function p is a disjoint sum of submodular functions A* for zero 
singular nodes i (see (4.2)). So the exchange capacity k{-]u,v) for a pair of u and v can 


take positive if u,v belong to Ui for some zero singular node i, and is equal to the exchange 
capacity for submodular function A* on a 6-element set Ui. Hence this can be computed in 
constant time. The number of nodes of Fp^r is at most 6n -|- 2, and the number of edges is at 
most 2m -|- 8n. Thus step 2 is done in 0(MSF(n,m, 1)) time. □ 


Proof of the main result (Theorem l.lD . As in Section construct nondegener¬ 
ate instance {N,S). Th en d(Fo) < 0(mlog fe). Apply the dual descent algorithm for 
{N,£). By Theorem 4.5, we obtain an optimal potential (p, r) and an optimal multiflow 
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/ in 0((m log/i:)MSF(n, m, 1)) time. As we have shown in Section]^ / is also a maxi¬ 
mum multiflow, and an optimal potential {p, r) for the original instance is obtained from 
(p,f) (in 0{nm\ogk) time). Then r is a half-integral optimal solution of LP-dual (1.2). In¬ 


deed, Yli£V\S c{i)r{i) is equal to the maximum flow-value. The feasibility of r follows from 
Eiey(P)\s = Eiiei?(P)(^(*) +^(i)) ^ J2ijeEiP)diPi,Pj) > 2 for every 5-path P. 


Proof of Lemma 4.3, Let X be the unique minimal minimum (a'’',o )-cut (of finite cut 
capacity). By Lemma 


2.2 


A is a transversal, implying (cl). Also X cannot meet {s^, for 
any terminal s. Otherwise, X contains s~ (and does not contain s^). However the deletion 
of s~ from X does not increase the cut capacity, contradicting the minimality. 

Consider a zero flat node i. Suppose that G X. If 0 A {k' ^ k), then the deletion 
of from A decreases the capacity, contradicting the minimality. Thus we have (c3). 

Consider a zero singular node i. Suppose that A n l/j / 0. If A n = {*^}, then A n Lj 
is of type 1, 3, or 5, the change A —)• A \ {i^} preserves the type at Ui, and the cut capacity; 
contradicting the minimality. Similarly, if A n C [/~, then the change A —)• A \ U~ 
preserves the cut capacity; a contradiction. Thus X r\Ui % Uf. So suppose that U~ n A = 0 
or U~ \{*^}. If U~ n A = U~ \{*^}, then necessarily = {i^}- Suppose that U~ CiX = 0. 

If |A n U^\ > 2 (type 1), then the change A —)• A U (Procedure (B)) does not increase 

the cut capacity (and keeps A being a transversal). In particular, the resulting A n I/j is one 
of the patterns in (c3). 

Consider a positive singular node i. First we show that |An {*)”) *3 }l ^ Suppose to 

the contrary that A n 17* = {i^}- The the change A —)• A \ {i)", fg } pi'eserves the capacity; a 
contradiction to the minimality. Also A n Uf = {ig } is impossible, and | A n {^i”) *2 > ^ 3 ^}l — ^ 
implies ig G A. Thus the pattern of A n U[~ is one given in (c2). Next consider A n . If 
A contains two nodes in ig"}, and necessarily A n U~ = 0; the change A —)■ A U 

(Procedure (B)) keeps A being a transversal, and does not increase the cut capacity. If A 
contains then A contains at least two nodes in , reduced to the case above. Thus, after 
Procedure (B), the resulting cut is a normal minimum cut, as required. 


Proof of Proposition In P* MZi*, vertices in B are colored black, and vertices in W 

are colored white. Other vertices have no color. For (p, r) G F* Kl Z* and p' G P* with p ^ p', 
let {p,r)\^pi := {u,r — 1/2) for the unique neighbor u oi p with d{p,p') = d{u,p') + 1/2. In 
addition, if {p, r) ^ B U W, then let (p, r)^pi := (u, r — 1/2) for the (unique) neighbor u of p 
with d{p,p') = d{u,p') — 1 / 2 . 

Lemma 4.6. For (p, r), {p',r') G F* lEI Z* with d{p,p') > 1, d{p,p') — r — r' is even if and 
only if one of following pairs has the same color: 

{{P,r),{p ,r')), ((p,r)\p/, (p',r')), ((p, r), (p, r')\p), {{px)\p',{p\r')\p). 

Proof. For (p, r), (p', r') G B U W, observe that d{p,p') — r — r' = d{p,po) -|- r -|- d{p',po) + r' 
mod 2. Hence d{p,p') — r — r' is even if and only if (p, r) and {p\ r') have the same color. Also 
observe that d{p,p') — r — r' does not change when (p, r) is replaced by (p, r)\p/. The claim 
follows from these facts. □ 

Lemma 4.7. There is no edge between Ujr and Ux- In particular, any normal cut X is 
decomposed into F-normal cut Xj := A n Bj and X-normal cut Aj := X nUx such that 

c{6X) - c(<5{a+}) = c(5A^) - c(5{a+}) + c{dXx) - c(J{a+}). 

Proof. Edge i^iy (or i'j^i^i) in Vp^r belongs to 17* for node i of colored {pi,ri), and hence 
belongs to Ujr or Ux. So consider an edge i^j^i for distinct i,j. Then d{pi,pj) > 1 (since 
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Qij > 2) and d{pi,pj) — ri — rj is even (since aij is even). By the previous lemma, if both 
{pi,ri) and {pj,rj) are colored, then they have the same color. Suppose that {pi,ri) has no 
color. If {pj,rj) has a color, then {pi,ri)\p. = {pi^*k-,'fi — 1/2) has the same color. If {pj,rj) 
has no color, then {pi,ri)\pj = {pi^*k,ri - 1/2) and {pj,rj)\p^ = (pj^*k>,rj - 1/2) have the 
same color. Consequently G Ujr or \^3k' ^ ^ 2 : for all cases. □ 

Let Fp^^ (resp. 2j))^) denote the subset of Fp^r (resp. Tp^r) consisting (p, r) with r* > 0 for 
i = 1,2,... ,n. Proposition |4.4| follows from the above Lemma |4.7| and the following. 

Lemma 4.8. (1) The map X — )■ {p,r)^ is a bijection between F^^j. (resp. and the set 

of all F-normal cuts (resp. X-normal cuts). 

(2) For an F-normal or X-normal cut X, it holds 

9n,£{{p, r)^) - gN,£ip, r) = c{dX) + p{X \ {a+}) - c((5{a+}). (4.4) 


Proof. (1). We claim that Xi := {X (1 Ui \ X : J^-normal cut} and Fp.^.,.. are in one-to- 
one correspondence by (4.3). Suppose that (pi,ri) G W, which implies Ui C Ux. Then 
^PiXi ~ ~ {(PiXi)} s-iid Xi = {0}. Thus the claim is true. Suppose that {pi,ri) G B, 

which implies Ui C Uj^. The claim can be seen from Figure]^ Suppose that (pj, r^) 0 B U IF; 
necessarily i is positive. We can assume that (pj_>.*i, r*-1-1/2) G B (and (pj^*i, r* —1/2) G IF). 
Then 7’+,., = {(pj,ri),(pj^n,r’i-l/2),(pi^.2,ri-hl/2)}. Since e Ux a.ndi^,if GUjr, 

it holds that T) = {0, {if}, {*^}}- Thus we have the claim, implying the statement (1). 

(2). Let {p',r') := {p,r)^. We first show that g{p',r') < 00 if and only if c{6X) is finite. 
Suppose that g{p',r') < 00. Pick if ^ X and edge iff^i- We show j^, G X. Recall that 

d{Pi,Pj) — ri — rj — aij = 0 and aij > 2. Hence d{pi,pj) > 2. Since if G X, the change 

{Pi, fi) —>■ (p(, r() increases d{pi,pj) — ri — rj by one. Necessarily the change {pj,rj) —)• {p'j,^^) 
must decrease d{pfpj) —r[ — rj by one. This means that {p'j,^)) = {PjUj + 1); (Pj^kiUj) or 
{pj^*k',rj + 1/2) must hold. For all cases, X contains jy. 

Suppose that c{6X) is finite. We show that g{p', r') < 00. Here r' > 0 is clear. We need to 
show that fi(p(,p' )—r'—r' —ajj < 0 for each edge zj G E. Suppose that d(pi,pj)—rj—rj — ajj = 
0. Namely ij G Ep^r and d{pi,pj) > 2. In this case, the argument is the same as the 
above. Indeed, suppose that the change (pi,rj) —>■ (p(,r() increases d{pi,pj) — ri — rj by one. 
Then X contains if with pj G F*, ^, and hence contains jy with pi G F*_ y. The change 
{pj,rj) -)■ {p'j,r'j) decreases dfpfpj) -r[- rj by one, implying dfpfp}) -r[- r'- - aij = 0. 

Thus it suffices to show that dfpfp'j) — r\ — r' — aij = 1 cannot occur. Otherwise, 

d{pi,Pj) -Vi- rj - aij = -1, d{pi,pj) -ri- rj is odd, and (p', r') ^ {pi, rf), {p'j,r}) ^ {pj,rj). 

and 
has 


If both {pi,ri) and {pj,ri) have colors, then these colors are different (by Lemma |4.6[ ) 
{p'i,r'i) = {pi,ri) or {Pj,'i"j) = {pj,rj) must hold; a contradiction. Suppose that {pi,ri 
no color. Then {pfr'f) = {pi,ri)^p. must hold. If (p},?’') has a color, then this color must 
equal the color of (pfr'f) (since {pi,ri)yp. and {pi,ri)\p- have different colors), and hence 
{p'jU'j) = {Pj,rj); a contradiction. If {p'jX'j) has no color, then {p'-,r}) = {Pj,fj)^pi must 
hold; but this is impossible since colors of {pi,ri)^p. and {Pj,rj)^p- are different. 

Finally we show the equation (4.4). The left hand side is equal to f2iev\s 
the right hand side is equal to 


{c{Ain6X)-2cii)}+ ciAinsx)+ 4*(^nt/i) 

iipositive zizerofiat zizero singular 


One can veri fy fr om the network construction (see Figure]^, definitions of {p,r)^ (see (4.3)) 
and Z\* (see (2.7)) that 2ci(r' — ri) is equal to c{Ai n dX) — 2c(z) if z is positive, c(4 n SX) if 
z is zero flat, and A}{X nUi) if z is zero singular. Thus we obtain the equation (|4.4[). □ 
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A Appendix 

For X G r* MZ* , the first and second components of x are denoted by Xp and Xr, respectively, 
i.e., X = {xp, Xr)- 


Proof of Theorem 3.5, First consider the case where F is a path (with no ends). Then 
r*MZ* is isomorphic to the product of two zigzagly-oriented paths. In this case, (F*IEIZ*)" is 
identified with the product (Z*)^"' of 2n paths, and L-convex functions on (E* Kl Z*)^ coincide 
with alternating L-convex functions in the sense of HZ! Also Doo is equal to d in HZI. Then 


Theorem 3.5 was shown in HU Theorem 2.6]. In particular, the following holds: 


Proposition A.l. For x G dom^r and a minimizer x' of g over fFx with g{x') < g{x), if 
g{x) = minygi^ g{y) or g{x) = g{y), then it holds 


Doc{x',opt{g)) = Dooix,opt{g)) - 1. 


(A.l) 


We show that this proposition holds for a general tree F. Pick an arbitrary z G opt (< 7 ) 
with T)oo(x,opt( 5 )) = Daoix,z). Since (x')p and {xi)p are equal or adjacent in F or F*, 
there is a path Pi in F* containing (xj)^, (x()p, and {zi)p. Then z,x, and x' are points in 
nr=i -Pj Kl Z* ~ (Z*)^*^, and the restriction g' of g to 0?=! -P* ^ (alternating) L-convex. 

Thus Proposition A.l is applicable to g',x,x', and hence Dao{x' ,opt{g')) = Doo{x,opi{g')) — 
1 = Daoix,opt{g)) — 1. By Dadx',opt{g)) < Doo{x',opt{g')), the equation (A.l) holds. 

Theorem [3.51 is proved as follows. By the description of the steepest descent algorithm, 
it holds that g{x^) = minj^gj^. g{y) or g{x^) = min^gj-^. g{y) for all i > 0. Thus m — 
1 = opt( 5 )) holds. Since Dao{x^,opt{g)) < Dao{x,opt{g)) + 1, we obtain m < 

Dooix,opt{g)) -b 2. In addition, if g{x) = min^gi^ g'(y) or g{x) = min^gyr^ ^(y), then m = 
Zloo(a ^5 opt(9')) holds. 
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Proof of Proposition |3.6| for a positive cost. For an even positive integer a(> 2), define 
h : {r** K Z** f ^ R by 


h{x,y) := d{xp,yp) - Xr - Vr - a {x,y e F* 


Lemma A.2. For x,y,x',y' G F*MZ*, it holds h{x,y)+ h{x',y') > 2h{{x + x')/2, {y + y')/2). 

Proof. A classical result |7l Lemma 3] in location theory says that the distance function on a 
tree is convex. Thus it holds d{p, q) + d{p\ q') > 2d{{p + p')/2, {q + q')/2) for p, q,p', q' e F*. 
The inequality immediately follows from this fact. □ 

Lemma A.3. For x,y £ F** Kl Z**, if h{x,y) < 0, then h{\x], [y]) < 0 and h{[x\, [yj) < 0. 

Proof. Let A := h{\x~\, [y]) — h{x,y) and A' := [yj) — h{x,y). Then it holds that 

A, A' G {—1, —1/2, 0,1/2,1}. Hence we can assume that h{x, y) G { — 1/2,0}. By a >2 and 
h{x,y) > —1, it must hold d{xp,yp) > 1. Then there are two edge-disjoint paths in F**, one 
containing \x'\p,Xp, [xjp and the other containing |'y]p,yp, Yy\p. Also they are contained by 
a single path. From this, we see d(|'x]p, \y'\p) — d{xp,yp) = d{xp,yp) — d{[x\p, [yjp). Also 
\x~\r — Xr = Xr — Yx\r and [ylr — yr = Ur — Yulr- Hence we obtain A' = — A. Notice that 
both h{\x'\, [y]) and [yJ) are integers. If Z\ = 1/2, then h{x,y) = —1/2, and hence 

H\x], \y]) = 0 and h{[x\,[y\) < 0. 

So it suffices to show that h{x,y) = 0 and Z\ = 1 cannot occur. Suppose not. Then 
both /i(|'x], [y]) and h{[x\, [yJ) are odd. If both [x] and [y] (or [xj and [yJ) are colored. 


then they have different colors (Lemma 4.6), x = [x] = [xj or y = [y] = [yJ holds, and 
A = 1 is impossible. Hence we can assume that Xp,yp G F** \ F*, and ([x], [yJ) £ W x B 
or (W, fyl) e B X W. Then d([x]p, [y]p) - d{xp,yp) = d{xp,yp) - d([xjp, [yJp) = 1/2, 
d{\x]p, [yjp) = d([xjp, [yjp) = d{Xp,yp), and \x'\r-Xr = [yjr -yr = Xr - [xj^ = yr - [yjr = 
— 1/4. Hence h([xj, [yJ) = /i([xj, [yJ) = h{x,y) = 0 (even). However ([xJ, [yJ) G IT x R or 
([xJ, [yJ) G R X VT implies that /i([xj, [yJ) or h{[x\, [yJ) is odd; a contradiction. □ 


We are ready to prove Proposition |3.6[ It is not difficult to see that the function on 
F* Kl Z* defined by x i—)• x^ if x^ > 0 and oo otherwise is L-convex. So it suffices to show 
that the function on {F* Kl Z*)^ defined by (x,y) i—)• 0 if h{x,y) < 0 and oo otherwise is 
L-convex. We show that h{x,y) < 0 and h{x',y') < 0 imply h{\{x + x')/2'], [(y-|-y')/2j) < 0 


and h{[{x + x')/2j, [(y -|- y')/2\) < 0. By Lemma A.2 we have h{{x + x')/2, (y -|- y')/2) < 0. 
By Lemma A.3[ we have the desired results. 
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